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Abstract 

Following the recent measurement of the acoustic peak by the BOOMERanG and 
MAXIMA experiments in the CMB anisotropy angular power spectrum, many anal- 
yses have found that the geometry of the Universe is very close to flat, but slightly 
closed models are favoured. In this paper we will briefly review how the CMB 
anisotropics depend on the curvature, explaining any assumptions we could make 
and showing that this skewness towards closed models can be easily explained by 
degeneracies in the cosmological parameters. While it is difficult to give independent 
constraints on the cosmological constant and/or different forms of dark energies, we 
will also show that combining CMB measurements with other observational data 
will introduce new and tighter constraints, like Qa > at high significance. 
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1 Introduction 



In the most general and simple inflationary scenario [1] the overall present 
energy density of the universe must be equal to the so-called critical energy 
density {Q = p/pc = 1, with pc = SHq/StiG). In fact, the condition for 
inflation {dinj > 0, and so d/dt{H~^f/ainf) < 0) is precisely that which drives 
f2(t) towards 1 in the Friedmann equation 

n{t) - 1 



a?{t)H^{t) 



during the inflationary period. This prediction, however, taken with the stan- 
dard CDM model of structure formation, is apparently in disagreement with a 
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combined set of observations, such as density- velocity galaxy field comparisons 
([14]), pairwise galaxy velocities ([15]), X-ray clusters temperature function 
evolution ([7], [8]) and velocity dispersions [^([9])], which point in favor of a 
low density universe [f] {^matter < !)■ In order to solve the discrepancy, open 
(Q = Qmatter < 1) inflationary models have been proposed ([2], [3] [4]). This 
type of model highlights limitations on the predictiveness of the inflationary 
scenario, which is supposed to have the advantage of removing any depen- 
dence on initial conditions from our present observable universe. However, 
even if this picture leads to a more complicated phenomenology, it generally 
determines f2 directly from parameters of the physical theory. 

On the other hand, another way to keep low-density models compatible with 
the simplest model of inflation is to introduce a cosmological constant. A, such 
that Q = ^matter + = 1- The prescucc of such a cosmological constant, 
which is compatible with all the above observations since they are practically 
insensitive to it, is also preferred by measurements of the magnitude-redshift 
relationship in high-redshift type la supernovae ([16], [17]). Nonetheless, this 
'natural' solution introduces the cosmological constant problem ([5]) that is 
perhaps even more acute in inflationary cosmology (see [6] for a review). Thus, 
an accurate determination of the present overall density parameter f2, even if 
not a panacea for the cosmological scenario, is at least extremely important in 
understanding which theoretical framework could explain the above conflicts. 

There is much experimental evidence for the presence of a peak in the CMB 
angular power spectrum ([18], [19]). Furthemore, with the recent release of the 
BOOMERanG-98 [20] and MAXIMA-1 [30] spectra, the shape and position 
of the peak has been detected with unprecedented accuracy. This result, apart 
from being in wonderful agreement with the standard scenario of primordial 
adiabatic fluctuations, has important consequences on the parameter ([33], 
[37], [38], [39], [40]). As previously noted and already well explained in the 
literature, the effect of the curvature is to change the relationship between the 
physical scales on the Last Scattering Surface (LSS) and the corresponding 
angular scales (see., e.g. [21], [22]). In an open universe, for example, the 
geodetics focalize in such a way that a particular angular scale will correspond 
to a greater physical scale on the LSS than the one expected in a flat model. 
The immediate result is a shift in the radiation angular power spectrum (the 
so-called C^'s), and thus a dependence of the position of the flrst peak ipeak on 
the curvature ([41], [23], [24], [25], [26]). The CMB power spectrum is therefore 
a powerful tool for the determination of the curvature and so of the overall 
energy density. 



A more conservative approach would say that the situation is rather unclear and 
the X-ray clusters data and theoretical modelling can be sufficiently large to prevent 
an unambigous exclusion of = 1 (see e.g. [10], [11], [12], [13]) 

From here on, with O we will indicate only the present value 
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Figure 1. The BOOMERanG and MAXIMA Q's data together with the most viable 
open, closed and flat models from present non-CMB observations. 

2 ft and the shift of the CMB Angular Power Spectrum 

In Fig.l, we show the recent data from the BOOMERanG ([19], [20]) and from 
the MAXIMA [30] experiments together with the predictions of the most vi- 
able open, flat and closed adiabatic models (from present non-CMB observa- 
tions) . It is quite evident that the open (closed) model predicts a first peak on 
smaller (larger) scales with respect to the flat model, which is in much better 
agreement with the data. 

In the literature, the dependence of ipeak on ft is often expressed as 



However, this approximation is not correct in A-dominated universes (see 
Fig.l, where the peak for the closed model is at ~ 130 instead of ~ 190) 
and does not take into account the further dependence of ipeak on other pa- 
rameters like the Hubble constant and the matter density Q,,„. So, in view of 
the recent ipeak = 197 ± 6 ([20]), further modiflcations to the above formula 
are needed (see also [27], [28]). Also, the width of the peaks and the inter 
peak distance vary as functions of ft (again see Fig.l) so a more complete 
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expression is needed to describe these effects. Tliis point is ratlier important 
because knowing the exact dependence of the CMB spectrum on fl will help 
us in understanding the shape of the probability distribution function for this 
parameter and, ultimately, how well it can be measured independently. 

In order to do this, it is convenient to introduce the 'shift' parameter R defined 
as 

ipeak — ^ieakl ^ (3) 



where flat indicates the flat, pure-CDM, A = model. The use of the R 
parameter is more appropriate than the conventional (,peak because it has a 
clearer geometrical dependence. 

As usual, let us assume the metric of space-time to be of a Friedman-Lemaitre- 
Robertson-Walker (FLRW) form with curvature k: 

ds^ = a(t)^[-dri^ + -fijdx'dx^] (4) 

with 

-iijdx'dx^ = dr^ + x^(r) {dO'^ + sin^^^^) (5) 



where the function x(r) depends on the curvature k and is r, sin(r) or sinh(r) 
for A; = (flat models). A; = 1 (closed models) and A" = — 1 (open models), 
respectively. The position of the first acoustic peak is determined primarly by 
the angle subtended by the acoustic horizon \ac at decoupling time, rjdec- The 
angle under which a given comoving scale A at conformal time rjaec is seen on 
the sky is given by ^(A) = X/x{Vo~Vdec)- As the harmonic number £ is inversely 
proportional to the angle 9, we have R — Oac/d^"'*, with 9ac — CsVdec/ xiv) ^ 
where = 1/^3(1 + 3il,b/AD,rad{^ + Zdec)) denotes the adiabatic sound speed 
of the baryon/photon plasma at decoupling. It is possible to show that [22]: 



Vdec = 



rod + ^m/{Zdec + 1) - \/^rad) (6) 



where z^ec ~ HOO is the redshift at decoupling and flk = 1 — flm — ^a- 
Furthermore, we have: 



dz 



y-VO- Vdec - ^ / [^^(^ ^ + + ZY + Qa]V2 
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For a flat, flm = 1, I^a = 0, universe we have 



^fc"* = CsVdec/ iVO - Vdec) = Cs{\j Vlrad + l/(2dec + 1) " \l^d) (8) 



We then find, keeping constant VL^h"^ (and so c^) and VLmh? (see next section) 
[29]: 



R 



2 la 



X(2/) V fir 



(9) 



which is a quantity just dependent on 17^ = 1 — — and VLr. 




Figure 2. R = constant lines in the f^A — plane. 

In Fig. 2 we plot the contours at i? = const in the — plane, together 
with one dashed line such that fi = 1. As we can see, open models have in 
general R < 1, while closed models have R > 1. This is the usual result 
that in an open model the peaks are shifted towards greater i values (smaller 
angular scales) with respect to the flat model case, while for closed models we 
have the opposite effect. It is worth noting that the CMB angular shift is not 
linearly related to the curvature: lines at Q = const are not parallel with the 
R = const contours, but have multiple intersections, especially in regions far 
from flatness (again, see Fig. 2). 

In Section 4, we will build a likelihood distribution function L for fl, using 
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a Bayesian approach to compare the current CMB data with the theoretical 
predictions and then marginalizing over the remaining cosmological parame- 
ters. The above result implies that the probability distribution function for 
Q (if we assume a flat prior distribution on the remaining parameters while 
marginalizing) will always be 'skewed' so that fl will never be measured at a 
level better than 10 - 20%[3. 



3 The geometrical degeneracy 

With the parameter R fixed, the structure and position of the spectrum 
is dependent on 2 physical scales: the equality scale and the sound horizon at 
decoupling scale. These quantities are completely defined once we choose the 
abundance of cold dark and baryonic matter in our model, by the parameters: 
'^cdm = ^cdmh^ and cjfe = VLi,h? . The CMB spectrum also depends on the 
characteristics of the primordial inflationary perturbations. Assuming that we 
have already selected the primordial power spectrum of our model, both in 
shape (tilted or 'blue') and in nature (adiabatic, isocurvature, hybrid), the 
structure of the CMB angular peaks is completely determined by R, ujh and 
ijJm = <jJb + ^cdm- cdm This result has an important consequence: if we let fib, 
flcdm and h assume any value but a fixed uOm and Ub, the lines at R = const 
in the Qm — plane correspond to sets of degenerate power spectra with an 
identical shape on subdegree angular scales {£ > 30) [29]. 

In Fig. 3 we draw a set of degenerate models with R, Um and Uf, set to those 
of the cosmological concordance model, = 0.35, = 0.65, h = 0.65, 
Ub = 0.021, Us = 1. It is clear from the degeneracy of the models that, once 
A is included, the peak position is not directly related to fl. It is also clear 
that it seems impossible to obtain any relevant and independent information 
from the CMB on fl^: the Integrated Sachs- Wolfe effect on large scales could 
in principle break the degeneracy but cosmic variance and a possible presence 
of a gravity wave background make this effect difficult to disentangle. 

Our main results are then the following: 

• Lines at i? = const in the — flm plane correspond to sets of degenerate 
Ce power spectra. 

• Given a flat fl\ ~ 0.65 model, a degenerate closed model can be found by 
decreasing Qa and h and increasing fib and flcdm- 

• R, ojm and Ub are the most meaningful CMB anisotropy observables. 

^ Another interesting point, is that for fim ~ the lines at i?=constant converge 
towards fl^ = 1, but we must warn the reader that in the de Sitter solution the 
notion of open, flat or closed Universe becomes ambigous [40]. 
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Figure 3. Concordance model and degenerate models 

• The CMB spectrum is a useful tool for the determination of Q only if we 
live in a flat universe {R = const are not parallel to = const when fl 
greatly differs from one). 

• Assuming that the concordance model describes our real universe we expect 
that a likelihood analysis for Q, using only the CMB power spectrum and 
without including any external information about h, Q,rn and f^A, will always 
be skewed towards closed models. 



4 From CMB to n 



Let us now describe the standard tools for extracting the cosmological pa- 
rameters from CMB anisotropy observations. Here we will analyze the recent 
BOOMERanG97 ([19]), BOOMERanG98 ([20]) and MAXIMA ([30]) results. 
The power spectra from these experiments were estimated in 7, 12 and 10 bins 
respectively, spanning the range 25 < i < 785. In each bin, the spectrum is 
assigned a flat shape, i{i + l)Ci/2n = Cb- Following [31] we use the offset 
lognormal approximation to the likelihood L. In particular we deflne: 

- 21nL = {D'^ - D%^)MBB'{Df, - D%^,), (10) 
D^ = HC^ + xb), (11) 
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Mbb' = (C|^ + xb)Fbb'{C''^> + XB'), 



(12) 



where Cg (Cg') is the theoretical (experimental) band power, is the offset 
correction and Fbb' is the Gaussian curvature of the likelihood matrix at the 
peak. Of course, will depend on the various parameters of our cosmological 
model {uJm,ns, ...), and so it will be the likelihood function L. In order to 
compute the likelihood for a given parameter only a we can either marginalize 
over all the remaining parameters, namely carry out the integral 

Lmarg{.OL) ^ j Pprior{oL,'^)L{a,'n)d'U. (13) 

— * 

where H is a vector containing all the remaining parameters and Pprior is the 
prior probability distribution, or we can maximize i.e. for a fixed a find the 
^max wich maximizes 

Lmaxip^^ Pprtor{p^i^max)L{oi.iT\.jjia^x) ■ (-^^) 



The two methods, in general, agree at a level of ~ 10%. The maximization 
method is usually based on a search algorithm through the second derivative 
of the likelihood matrix ([32]). In this approach the Q spectrum is computed 
on the way, without sampling the whole parameter space. A different approach 
is based on building a database of C^'s on a discretized grid of the parameter 
space. L{p) is then obtained by maximizing and/or integrating the likelihood 
computed on the grid ([33], [34]). Of course, producing a grid of models can 
be quite computationally expensive even with the new and fast boltzmann 
codes like CAMB [43] or CMBFAST [42]. But this problem can be drastically 
reduced using morphing [35] or interpolation [36] algorithms. 

The definition of the database is rather important because it defines the in- 
our analysis and, in general, it is better to have this prior as 
flat as possible for each given parameter. This brings us to the choice of the 
variables in which the database must be sampled. As we saw in the previous 
section, the CMB anisotropics are mainly sensitive to the physical variables cui 
and R, so sampling in those variables will avoid degeneracies. Furthemore, the 
physical baryon density oob = ^bh'^ is well determined by independent mea- 
surements like primordial nuclide abundances, so this is the optimal choice for 
extracting information about this parameter (without involving complicated 
Jacobian transformations) or assume external priors for it. For the same rea- 
son, extracting confidence limits on parameters like h can be a little more 
elaborate with this sampling, being the database in ~ /i^. 

Another possibility is to sample the database in cosmological variables, like 
and h. Of course, this will introduce degenerate models in our database 
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but this sampling has the advantage of obtaining direct constraints on the 
commonly used parameters and with flat prior distributions. In most of the 
recent papers, either a 'hybrid' variables approach, with sampling in Ui and 
1—flm — ^A, or the cosmo/o5'«ca/ variables approach has been used. Here we will 
choose the database approach, sampling the parameter space in cosmological 
variables as follows: Qm = 0.1, ...,1.1; = 0.0, ...,1.0; Qb = 0.01, 0.25; 
h = 0.25, 0.95 and Us = 0.50, 1.50. We will not consider the possibility 
of high redshift reionization of the intergalactic medium Tc >> 0, a gravity 
waves contribution or the effect of massive neutrinos. 



5 Removing the geometrical degeneracy: Results 



B98+B97+MAXIMA-1+COBE 

Full database, age > lOGvr prior 




68% c.l. 



: 0.23 <£1^< 0.75: 
:0.26<£1^<0.83 
iO.98 <£!< 1.17 



Figure 4. Likelihood functions with age prior. 

In Fig. 4 we plot our likelihood contours for Q, Qm and using just the 
intrinsic internal priors of the database plus the quite reasonable age prior 
tuniverse > lOGj/r. As wc cau scc, the likelihood for fl is skewed towards closed 
models but is consistent with flatness. It is rather important to note that this 
skewness is largely due to the i?-degeneracy which makes 'more' closed models 
compatible with the observations. The likelihood for Qa and flm are quite flat 
due to the geometrical degeneracy but they nonetheless feel border effects from 
the database priors. The likelihood for fl starts to be in even more agreement 
with flatness when a Gaussian prior h = 0.65 ± 0.2 is assumed as in Fig. 5. 
This clearly shows that most of the degeneracies in the h < 0.5 region are well 
removed by the prior. Including a prior flm = 0.3 ± 0.2 (Fig. 6) as suggested 
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B98+B97+MAXIMA-1+COBE 
Full database, age > lOGyr prior, 0.45<h<0.85 
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Figure 5. Likelihood functions with h = 0.65 ± 0.2 and age priors. 



B98+B97+MAXIMA-1+COBE 

Age prior, 0.45<h<0.85. 0.1 < n ,, < 0.5 
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Figure 6. Likelihood functions with = 0.3 it 0.2 prior. 

by the majority of measurements, shrinks the likelihood towards = 1 and 
gives a strong determination for the cosmological constant Qa > at Aa level. 
The complementarity with the supernovae type la measurements is even more 
clear in Fig. 7, where a combined CMB+SnIA analysis gives > at more 
than 7a and VL = 1 with a few percent uncertainty. 
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B98+B97+MAXIMA-1+COBE+SN 
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Figure 7. Likelihood functions for a combined analysis with Supernovae. 



6 Conclusions 



The BOOMERanG and MAXIMA data support the main prediction of the 
inflationary paradigm: that the geometry of the universe is flat. The small 
deviations towards closed models reported in various analyses ([37], [39]) can 
be easily explained by the degeneracies in the cosmological parameters, which 
make more closed models compatible with the data. These conclusions are 
considerably strengthened by the inclusion of other cosmological data such as 
measurements of the Hubble constant, the overall matter density Qm and the 
accelerating expansion rate indicated by observations of distant Supernovae. 
At the same time, fl\ and other forms of 'dark energy' cannot be well de- 
termined by the CMB data alone, in spite of their high precision. This does 
not mean that CMB measurement are not useful in the determination of such 
parameters: combining the CMB data with constraints from observations of 
large-scale-structure and from observations of SN-Ia increases the extent to 
which A can be quantifled, with f^A > at ~ 7a. Future data from the 
PLANCK and/or SNAP satellites will hopefully enable us to resolve the 'dark 
energy' puzzle. 
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